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Abstract 



^ (— i Pseudo-Hermitian operators appear in the solution of Maxwell's equations for stationary 

non-dispersive media with arbitrary (space-dependent) permittivity and permeability tensors. 
We offer an extension of the results in this direction to certain stationary dispersive media. 
In particular, we use the WKB approximation to derive an explicit expression for the planar 



time-harmonic solutions of Maxwell's equations in an inhomogeneous dispersive medium and 
study the combined affect of inhomogeneity and dispersion. 
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5 ' 1 Introduction 



A linear operator A that acts in a Hilbert space is called pseudo-Hermitian, if one can find an 
invertible Hermitian operator t] satisfying = rjHr] -1 , [T]. This notion of pseudo-Hermiticity 
arises naturally in the study of non-Hermitian Hamiltonian operators, such as H = — + ix 3 , 
that admit a real spectrum [2]. It has also interesting applications in many different areas [3]. 
Perhaps one of the most remarkable of these is the appearance of pseudo-Hermitian operators in 
classical electrodynamics [3]. 

Consider the propagation of electromagnetic waves in a source-free medium with arbitrary 
(possibly space-dependent) permittivity and permeability tensors £ and fl. Then, the Maxwell 
equations read 

V-/3 = 0, V-5 = 0, (1) 
B + DE = 0, D-DH = 0, (2) 

where E and B are the electric and magnetic fields, D = £ E, H = fl B, an over-dot means a 
time-derivative, and D is the curl operator, e.g., DE := V x E. 
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For the cases that the medium is stationary, i.e., £ and ]l do not depend on time, we can reduce 
Eqs. (ED to 

B(r,t) = B (f) - [ DE(x,r)dr, (3) 
J o 

£ + fi 2 £ = 0, (4) 

where B (r) := -B(r, 0) and 

Q 2 := £ D/i D. (5) 
In view of (TjQ), the initial- value problem for Maxwell's equations admits the formal solution [3]: 

E(r, t) = cos(nt)E (f) + IT 1 sin(fit)£ (r), (6) 

where £ (0 := £(f, 0), and E'o(f) := £'(r, 0) = £ D /i S (r). 

f> O O — 1 <r^ — 1 

We can view D, £, fl, £ and /i as linear operators acting in the vector space of complex- 
valued vector fields F : R 3 — > C 3 . Endowing this space with the L 2 -inner product, -< F,G >-: = 
J r3 F(r)* ■ G(r) dx 3 , we find a Hilbert space % in which the curl operator 2) acts as a Hermitian 
operator, i.e., for all F,G belonging to the domain of 2), we have -< F,DG >-=-< DF,G y. 
For a lossless (and gainless) medium, the permittivity and permeability tensors (and their inverses) 

also define Hermitian operators acting in H. This in turn implies f2 2 ^ = £Q, 2 £ , i.e., as an 
operator acting in "H, Q 2 is a pseudo- Hermitian operator. 

In fact, because for a lossless medium £{r) is a positive-definite matrix for all r, £ is a positive- 
definite operator, and fl 2 belongs to a special class of pseudo-Hermitian operators, called quasi- 
Hermitian [5]. A basic property of quasi- Hermitian operators is that they are related to Hermitian 
operators via similarity transformations [B]. For example, we can relate fl 2 to a Hermitian operator 
h according to 

n 2 ^h-.= ¥ q 2 T l = r 1 © T\ (7) 

This observation has been used in [1] to devise a method of solving the initial-value problem for 
Maxwell's equations. This method applies to arbitrary stationary inhomogeneous and anisotropic 
media, but it ignores the effects of dispersion. The purpose of the present article is to extend the 
approach of [3] to stationary dispersive media. 



2 Time- Harmonic Solutions and Dispersion 

In standard textbook discussions of the propagation of the electromagnetic waves in dispersive 
media, one usually begins by assuming a harmonic time-dependence for the electromagnetic waves, 

E(f, t) = e- iujt S{f) , B(f, t) = e~ iuJt £{f) , (8) 

and considers the possibility that the permittivity and permeability tensors £ and /l depend also 
on the frequency u. In Eqs. (jHJ), $ and 0B are time- independent complex- valued vector fields that 
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in view of ([T]) are subject to the constraints 

V ■(£#) = (), (9) 
V-J = 0. (10) 

Because according to ([3]) we can express the magnetic field in terms of the electric fields, we will 
confine our attention to the dynamics of the latter. 

First, we consider the case of non-dispersive material. 

Inserting the first of Eqs. (JSJ) in (jlj) yields 

Q 2 i=u 2 i. (11) 

This reveals the interesting fact that the solutions of Maxwell's equation having harmonic time- 
dependence are actually obtained from the eigenf unctions of the pseudo-Hermitian operator Vt 2 . 

In general, the eigenvalues of Q 2 are highly degenerate, and the imposition of (Q does not lift 
this degeneracy completely. In the case that the wave propagates in vacuum, we have £ = £ 1 , 
/i = Vol, and fl 2 = c 2 ® 2 . Therefore, the eigenvalue equation (ITT]) subject to the condition (J5J) 
coincides with that of the Laplacian subject to the same condition, 

- V 2 / = k 2 ^ , V • S z = 0. (12) 



Here k := oo/c and the superscript refers to the vacuum solution. As is well-known, Eqs. (TT2 
admit the following complete set of plane-wave solutions 



(2^p7 2 ~* 



<4V(0:= 7^575 Sf> (13) 



where k and et are any pair of orthogonal unit vectors, e- := k x e~~ , and k := kk. The unit vectors 
k and et play the role of the degeneracy labels for the eigenvectors. Substituting (TT3"j) for ${f) in 
(jSJ), we find the following solutions of Maxwell's equations having harmonic time-dependence 

p i(k-r— u)t) 

The general solution of Maxwell's equations in vacuum is a superposition of these plane wave 
solutions, 

' c + {k\ ^)El &t(vK) (x, t) + c.{k\ <Pi)^ Vs) (iB, t)\ . (15) 



p PZTY 

E (x,t) = d 3 k 

Jr 3 Jo 



Here for each k ^ we fix a polar coordinate system (p^ip£) in the plane normal to k, ^{f) is 
the unit vector specified by the polar angle fj: in this plane, er(^) = k x ef (<£]*), and c± are 
complex-valued mode functions. 

Now, consider the more general case that the medium is locally inhomogeneous and/or anisotropic. 
In this case, we will suppose that the degeneracy structure of the eigenvectors of Q 2 is the same as in 
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the case of vacuum. Therefore, we denote these eigenvectors by ~± , and express the corresponding 

' k 

solutions of Maxwell's equations with harmonic time-dependence (JS} as 



Ets^t) = e-^ k Mr). (16) 



-itut 

~k "'~k 



It is usually customary to take the general solution of Maxwell's equations to be a superposition of 
these time-harmonic solutions. 

For the case that the medium is lossless, Q 2 is a quasi-Hermitian operator with a complete 
set of eigenfunctions $tg±- This ensures the existence of an eigenfunction expansion for Q 2 and 

' k 

provides the mathematical justification for the claim that the most general solution E of Maxwell's 
equations is a superpositions of E% & ±, 



- j 



E(f, t) = I d 3 k [ d n \c + (k, y % )E m{ } (x , t) + c.(k, n )E %A - { Ax, t)] . (17) 
Jr 3 Jo k ' k j 

Next, consider the case of a stationary dispersive medium. Then, one usually accounts for 
the cu-dependence of £ and fl in determining the time-Harmonic solutions E^.± and considers 

' k 

the superposition of these solutions to describe waves having more general initial conditions, e.g., 
localized wave packets. This requires solving ( TlT]) . Note, however, that the operator Q 2 depends 
on 00. Therefore, this equation can no longer be viewed as a standard eigenvalue equation. Yet one 
can still perform the transformation (J7J) to determine an w-dependent operator h and express the 
solutions of Eq. (TTTT) in terms of those of 



hip = u 2 ip, 

according to 



S = £ >. (19) 

In the absence of dispersion, h is a genuine Hermitian operator with a complete orthonormal set 
of eigenvectors. Ref. [1] uses this observation to evaluate the right-hand side of the formal solution 
of Maxwell's equation given by (E])- Here we explore the consequences of allowing £, /I, and h to 
depend on oj. Because a direct extension of the results of [1] to general dispersive media is beyond 
our knowledge, we will confine our attention to separable cases where 



£(r,u) = e 1 (u)£ 2 (r), /i(f, u) = ^(u) /J, 2 (r) , (20) 

Ei, ui are positive real- valued functions of to, and £2, fl 2 are positive 3x3 matrix- valued functions 
of f\j The separability condition (120]) is not quite realistic, but as we will see in the following, it 
provides us with a concrete analytically treatable class of toy models. Furthermore, we can use it 
to establish the completeness of the eigenfunctions of h. 



1 In terms of the transformed permeability and permittivity tensors, £ (t,t' ,f,f') and fl{t,t' ,r,r'), 

<-> 

that appear in the constitutive relations [7]: D(r,t) = J_ dt' J H3 d 3 f' £ (t,t' ,r,r')E(r' ,t'), 
B(f,t) = fl^dt' J R3 d 3 r' jl(t,t',r,f')H(r',t'), The separability condition (J20J) takes the form £{t,t',r,r') = 
£i(t - V)£i{f)6{?' - r), jl(t,t',r, r') = fat - t')fa?)S(r' - r). 
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For the cases for which (j5DJ holds, the u- dependence of f2 2 and h factors, and we have 

a 2 = m(u)- 2 nl, h = ni (uj)- 2 h 2 , (21) 



where 

11 1 i i 1 i 

nx(u) : = v'eifwKH, ftj := £ 2 D fl 2 D , h 2 := £ 2 ilj £ 2 2 = £ 2 2 D jj, 2 D £ 2 2 . (22) 

As seen from (|22|) . ft,2 is a Hermitian (w-independent) operator acting in "H, and the eigenfunctions 
ip of h with eigenvalue w 2 coincide with the eigenfunctions ip 2 of h 2 with eigenvalue a; 2 n 1 (u;) 2 . 
Because h 2 is Hermitian, the latter form a complete ortho normal basis of "H. This in turn implies 
that the eigenfunctions of h will also form a complete orthonormal basis of % provided that the 
function 

f(u) :=wni(w) (23) 

is a monotonically increasing function of u; and we normalize the eigenfunctions properly. This 
together with ffl9|) constitute a mathematical basis for expanding the general solution of Maxwell's 
equations in terms of time-Harmonic solutions for dispersive media satisfying (1201) . 



3 Planar Waves in an Inhomogeneous Dispersive Medium 

Consider an isotropic but inhomogeneous dispersive medium defined by 

£ = e(z, w) 1 , /i = fi(z, (24) 

where e and fi are positive real-valued functions. Suppose that a linearly polarized planar wave 
propagates along the z-axis in this medium in such a way that both the electric and magnetic fields 
are independent of the x- and ^/-coordinates. The initial data have the form 

E (r) = E (z) = £ (z)i, S„(r) = B (z) = B (z)j, (25) 

where i and j are respectively the unit vectors along the x- and y-axes. 
Introducing p := —i4-, we can express the operators VL 2 and h as 

Vt 2 = e(z, w)~ 1 p/i(z, oj^p, h = e(z,<jj)~ l l 2 p ii(z,ijj)~ l p e(z,ijjy l l 2 . (26) 

Therefore, similarly to the non- dispersive case considered in [1], Eq. (1TB]) reads as 

1 d 

y/e(z,u) dz 

Here we have dropped the vector sign from if), because it is always parallel to the x-axis; ip(z) = 
if)(z)i. 

Eq. (I27p is the time-independent Schrodinger equation for a point particle moving along the 
z-axis and having a position- and energy-dependent-mass. Following jl], we will employ WKB 
approximation to solve this equation. This yields solutions of the form 



d 



li{z,uS) dz \y/ £ (z,u) 



cu 2 ijj{z) 



(27) 



where v(oS) are nonzero complex normalization constants, and 



u(z, lo) :- 



v(z, oj) 



y/£(z,U))(Jt(z,U)) 



(29) 



Substituting ( 128]) for ip in (|T9!) and using (124}) . we find a set of solutions S u of Eq. (ITTi) . Using these 
in OH]) yields the time-harmonic solutions for the system. In order to address the completeness of 
these solutions, we should address the completeness of (I28j) . We are able to do this for the separable 
cases where 

e(z,u) = e 1 (u)e 2 {z), fJ>{z,uj) = ^{^^(z). (30) 

The result is the following WKB-approximate eigenfunctions of h. 



D if(u))u 2 (z) 



2irv 2 (z) 

where / is the function given in (|23|) . a prime stands for a derivative, and 



v 2 {z) :-- 



^e 2 {z)n 2 {zy 



u 2 [z) :-- 



Ms)' 



(31) 



(32) 



Under the assumption that / is a monotonically increasing function, we can easily establish the 
following ortho normality and completeness relations for ip^. 



ipu(z)*il)u,>(z) dz = 5(uj - oj'), 



'4 } Lu(z)'ip L o(z'y du = 8(z - z'). 



(33) 



Having obtained a complete set of eigenfunctions of h, we can construct time-Harmonic WKB- 
approximate solutions of Maxwell's equations. In view of (jHJ) and (1191) . these have the form 



E w {z,t) 



e~ wt tl) u {z)i 



y/ei(u)e 2 {z) \e!(u)e 2 (z) 



fj.i(u)(i 2 {z) 



1/4 



1 

2^ 



1 + 



exp {iui[ni{uj)u 2 {z) — t]} i. 



(34) 

Note that the condition of the validity of the WKB approximation is the same as in the non- 
dispersive case (ei = n\ — 1) discussed in [1], namely 



2v 2 v% - v' 2 2 2fi2f4 ~ 3/i' 2 2 



2vl 



2(4 



(35) 



where we have suppressed the ^-dependence of v 2 and \i 2 for simplicity. 

In light of (I35p . the general WKB solution of Maxwell's equations with initial conditions of the 
form ( )25|) reads 



E{z,t) 



c(cu)E u (z, t) dcu, 



(36) 



where the modulus of the mode function c(u) is supposed to be negligibly small for all u violating 
(135]) . We can use f]33]) . ( 134"]) . and ( 136]) to express the mode function c(u) in terms of the initial 
electric field. This yields 



(u) = y/ei(u) j \/e 2 {z) ip u (z)*S (z) dz, 



(37) 



where Sq(z) := E(z, 0) • %. Therefore, the condition of the validity of the WKB approximation 
restricts the choice of the initial condition. 

If the medium happens to be homogeneous, i.e., £2 and /12 are constants, the WKB approxima- 
tion is exact, V2 is constant, U2 takes the form 7/2(2) = z/v2, and the time-harmonic solution (jBl 
coincides with a plane wave, Ae t ^ kz ~ ult \ where 




A : = 

k := ton\{uj)/v2 — us/s{uS)ijl{u). 
Therefore, if we consider the time-evolution of an initial plane wave: 

£n(z) = A e ik \ 



(38) 
(39) 

(40) 



we will not be able to observe the effect of dispersion. The situation is the opposite in an inhomo- 
geneous dispersive medium. Because of the z-dependence of v 2 and £ 2 , inserting (14"0~]) in (13~7]) gives 
a mode function c(oj) that is not proportional to a Dirac delta function. Consequently, the shape 
of the field changes in time; in an inhomogeneous medium an initial plane wave feels the effect of 
dispersion. 

For example, consider a nonmagnetic material (//i(w) = 1, ft 2 (z) = //o) with a Lorentzian 
inhomogeneity: 



e 2 lz 



£o[ 1 



a, 7 G 



(41) 



and an arbitrary dispersion relation given by £1(0;). Then to the first order in the strength of the 
inhomogeneity a, the mode function (|37|) corresponding to the plane wave (140]) takes the form 



where 



c(u) = 2t(w) <^ 5(k - &(u)) + 



2l(w) := A-t 



-7|fc-^)l( 3 fc-^(a;)) 
8(fc - ^(w)) 



1 Q 



'27T£o£lM 3/2 



1 



:= 



-0(a 2 ), 



(42) 



(43) 



2ei(w) 

The term in the square bracket in (142 j) reflects the dispersion of the plane wave (T40l) due to the 
inhomogeneity. It is interesting to see that this term is exponentially suppressed for the wave 
numbers violating the dispersion relation k = 



4 Conclusion 

In this paper we have examined the role of pseudo-Hermitian operators in the description of elec- 
tromagnetic waves propagating in a stationary, possibly inhomogeneous or anisotropic dispersive 
medium. 

In the absence of dispersion the properties of the pseudo-Hermitian operator Q 2 associated 
with the Maxwell equations lead to a powerful spectral method for solving the initial-value prob- 
lem for these equations pE]. In particular, for an effectively one-dimensional model, it yields an 
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explicit expression for the propagating wave provided that one can employ the WKB approxima- 
tion. When the medium is dispersive the same approach can not be pursued. Nevertheless, the 
pseudo-Hermitian operator Q 2 still generates the time-harmonic solutions as its eigenfunctions. 
This raises the question of the completeness of these eigenfunctions. In this article we addressed 
this question for the cases that the frequency- and the space-dependence of the permittivity and 
permeability tensors are separable. Although this condition was rather unrealistic, it allowed for 
a concrete implementation of our general method. In particular, for an effectively one-dimensional 
(planar) system with this separability property, we were able to construct a complete set of WKB- 
approximate eigenfunctions of Q 2 and study the combined effect of dispersion and inhomogeneity 
showing that unlike for a homogenous medium in an inhomogeneous medium a plane wave also 
undergoes dispersion. 
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